Josephson effect for superconductors lacking center of inversion 
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Abstract 

Due to the absence of a center of inversion in some superconducting compounds, a p-wave 
admixture to the dominant d-wave order parameter must exist. If time-reversal is also violated, 
an allowed invariant is the product of the d-wave, p-wave and an appropriately directed current. 



5-H 

^ ■ 

We show that this leads to an anomalous Josephson current for tunneling along the direction 

parallel to the axis of the p-wave component, where the current is the Meissner current and the 

Josephson loop current along the surface of the tunnel barrier. These ideas are applied to the 

the normal phase some experiments have detected a time-reversal and inversion symmetry broken 

phase. The effect is relevant also to heavy-fermion compounds which lack center of inversion due 
> 

to crystalline symmetry. 



PACS numbers: 



A microscopic state violating time-reversal and inversion symmetry and reflection sym- 
metry on some planes has been proposed for the pseudogap state of the cuprates |l|. The 
microscopic ground state associated with such symmetry breaking has spontaneous current- 
loops in the O- Cu-0 plaquettes in each cell from zero temperature all the way to the 
pseudogap temperature T*. Two classes of experiments ( ARPES with circularly polarized 
photons J2,l2| and polarized neutron diffraction j^|) on two different cuprate compounds have 
observed such a state below T*. These broken symmetries are expected to continue into 
the superconducting state, as indeed found in the experiments % |4||. This should then also 
affect the symmetry of superconductivity a,la|. Some heavy-fermion superconductors lack 
center of inversion in their crystal structure and similar superconducting states are expected 
for them p], [s| 

Quite generally, if time-reversal and inversion symmetry are broken, a term in the free- 
energy of the form 

ie(ipdi>px ~ c - c (!) 

is allowed, ipd is the dominant "d-wave" superconducting order parameter, ip px is a "p-wave" 
order parameter and the coefficient e is proportional to the order parameter associated to 
the symmetry breaking in the normal phase (for example the crystal structure lacking center 
of inversion p| or the spontaneous current-loops |l| in cuprates). x corresponds to the (110) 
direction or the (110) direction depending on the domain. We have introduced i explicitly in 
the coefficient so that e is real. In this case a superconducting state violating time-reversal 
and inversion symmetry and consistent with the proposed normal state is of the generic 
form ip{d xy ) + i5oip(p Xt y) in zero magnetic field, where the (110) direction in the crystal 
has been taken as the x direction. 5q = e/2ad, where the leading term in the free-energy 
for the d— wave component is al^l 2 . Such states have actually been investigated in con- 
text of some other materials whose crystal structure lacks inversion symmetry |2lla|- Some 
experiments have already been proposed to discover evidence for such a superconducting 
state in the Cuprates [5(. Here we will show that the Josephson effect between such super- 
conductors and a conventional s— wave superconductor with a particular orientation of the 
tunnel barrier with respect to the reflection symmetries of the superconducting state has a 
distinct signature from conventional Josephson effect. The clear observation of such effects 
will further substantiate the nature of the pseudogap phase for the cuprates and the nature 
of superconductivity in the relevant heavy-fermion compounds. 



Consider a surface of the superconductor perpendicular to the x-direction (110). Then 
for the state ip(d xy ) + iSoip(p x ), spontaneous super-current flows in the y-direction, i.e. along 
the surface, which is allowed. The alternative state ip(d xy ) + i5oip(p y ) has current normal to 
the surface and is not allowed. So the presence of a surface whose normal is the x-direction 
will naturally favor the p x admixture. Let us now consider that on this surface a junction 
is fabricated which allows a Josephson coupling to a conventional s-wave superconductor. 
We now wish to consider the modification of the superconducting state due to an applied 
magnetic field in a direction in the plane of the junction. This modification must come 
about because beside the invariant of Eq.(l), there also exists an invariant proportional 
to ijyip(p x )ip(d X y). This is because j y , the current in the ^-direction, is odd under time- 
inversion (as is iip(p x )) and the product of the three terms satisfies all inversion and reflection 
invariances. Such a term can be derived from the Ginzburg-Landau free-energy for the 
problem as follows: 

Let the magnetic field B be in the Indirection and adopt the London gauge in which the 
vector potential is 

A = Bxy. (2) 

The Free-energy density may then be written as 

F(^(4J,V'(Px),A) = a d |V'(4,)| 2 + ^|^(4,)r + /5d|V'(4 y )| 4 + a p |V'(Px)r+ (3) 

i€{il)(d xy )*ip{p x ) - ip{d xy )il){p x )*) + ie'{ip{d xy )*D y ip{p x ) - il)(d xy )D*ip{p x )*), (4) 

— * — * 

where D = — i(V — 2eiA/hc). This is simply the modification of the free-energy written 
down by Agterberg and Kaur jy| for the present geometry and in the presence of a magnetic 
field, e has already been discussed; e' has a magnitude of the order of the fermi-velocity 
and is independent of the underlying inversion/time-reversal breaking order parameter. On 
minimizing with respect to both the magnitudes and phases of d xy and p x , one obtains that 
the order parameter is 

ip = ip(d xy ) + i5tp(p x ), (5) 

where 

2e f x 

^(d xy ) = \ip{d xy )\ exp (i0 0x + i(f> B (x)), (f> B (x) = —y I B(x')dx' (6) 

and 

S=[6 - Oy]A d expi0 B (z)- (7) 



Here A^ is the magnitude of the major order parameter, 9q x = xee' jK R specifies the spiral 

n 

state deduced by Agterberg and Kaur p, 5 = (e(l + e' 2 / K R )) is magnitude of the bulk in- 
duced p— wave order parameter and the coefficient K R = K — e' 2 /a p and ( = e' '/ ' {2a p /S. 2 d K R ) . 
The length associated with the pitch of the spiral can be estimated to be much larger than 
the London penetration depth A. The spiral then produces negligible effects in Josephson 
tunneling and will be ignored by putting 9 = 0. 

The second part of Eq.(7) is equivalent to having the aforementioned term in the free- 
energy proportional to ijyipiPx)"^^^). There are two contributions to j y . One is just the 
Meissner screening current, which is uniform on the surface, 

— cB 
4 1J = 4^-exp(-z/A). (8) 

The other arises when a Josephson junction is constructed of this superconductor with 
another superconductor (of, for example the s-wave kind) and a magnetic field is applied 
in the junction as above. Then, the Josephson current density across the tunnel junction, 
J x (y) is periodic in y, 

■h (y) = Jxo sin(Bwy/$ + 7o) • (9) 

Here w is the effective thickness of the barrier in the x— direction and 70 is the phase 
difference across the barrier. It then follows, by continuity or by considering the gradient of 
the phase on the surface of either superconductor, that there exists also a periodic current 

3 ( y\y) = JM- (10) 

Such a periodic surface current is of-course present in every Josephson junction. Usually, it 
has no consequence. In the case of the d + ip superconductor, it changes the magnitude of 
the pa,— wave component periodically in the y— direction. Note that since this modulation is 
itself proportional to the Josephson current, it will have an effect on the Josephson current 
which is second order in the tunneling. We mention it for completeness although it has a 
negligible effect on our results and is dropped below. 

Collecting the results above, the wave-function of the d + ip superconductor near the 
surface of the Josephson junction up to a depth about A in x-direction is 

V(y) = Mfd(0) + iS y cos(0)) exp {iMy)) (11) 

= |A(j/)|expi(^ B (j/) + 0(y)), 



where 



A = A d yJf*(6) + 6*(y)co8*(6) (12) 

<fi = axctan(5 cos 8/ fd(6)) (13) 

5 = {6 + S 1 B), S 1 = Cc/(4vrA). (14) 

We may now calculate the Josephson current between the d + ip superconductor oriented 
along the x or (110) direction and a conventional s-wave superconductor. 

c-k/2 pL/2 

d6F(6) \ 

t/2 J-L/2 

F{9) gives the tunneling cone and L is the length of the barrier in the y— direction, i.e. along 
(110) . We can evaluate the integral to get 

-L/2 i-L/2 

dy (sin( 7o + (2irB\ L /$ )y)) + A 2 

-L/2 J-L/2 

where A x = A ^J 2 dOF(0)f d (6) and A 2 = A ^J 2 d6F(6) cos(0). We obtain 



J, ex/ dOF{6) [ dyAsm^o + ^-BXy + ^y)}. (15) 

J-n/2 J-L/2 n 



fL/2 f-L/2 

A x / dy (sin( 7o + (2irB\ L /$ )y)) + A 2 dy (5 cos( 7 o + (2nB\ L /$ )y)) , (16) 

J-L/2 J-L/2 



Li =Al sill(7o) || sill( ^) + A, (ft + « lB )cos( 70 )^sm(^)) (17) 

$ = BXiL is the total flux through the junction. This expression should be maximized with 
respect to 70 to get the observable Josephson current. The Ij thus obtained is 



h 

L 






(18) 



The first term in (|18jl is the usual term; It would be zero in a tetragonal crystal where 
fd(0) = sin(2#). The second term has a part due to the p-wave admixture; the more 
interesting part of it is proportional to B arising from the Meissner current, which cancels 
the flux $ in the denominator to give a part which as a function of B is oscillating with a 
constant amplitude. This last feature is remarkable since it will lead to a non decreasing 
amplitude for the Josephson currents with magnetic field for large fields. This happens due 
to quite different physics in squids and at a quite different period. 

The results of the evaluation of Eq. (18) are represented on figure 1, for a conventional 
d or s or d + s material (red dotted line), and a d + ip material (blue solid line). The main 
observation is that the critical current does not decrease at large fields but oscillates. For 



13 

C 

CD 

i— 

s— 
13 
O 

c 
o 
w 

-C 

Q. 
CD 
</) 
O 



03 
E 
X 



0.14 - 



0.12 - 



0.10 - 



0.08 - 



0.06 - 



0.04 - 



0.02 



0.00 



i — i — i — i — i — i — i — i — i — i — i — i — i — r 



l — i — i — i — i — i — i — i — i — i — i — i — i — r 



-15 



superconductor with center of inversion 
d+ip superconductor (one domain) 




-10 



-5 5 

Magnetic flux in units of 4> 



10 



15 



FIG. 1: (Color on line.) Josephson current patterns as a function of magnetic field. Red dotted line: 
Conventional Josephson pattern expected for a superconductor. Blue solid line: Josephson pattern 
expected for a single domain d + ip superconductor. The parameters used for the calculation in 
reference to Eq.(17) are A\ = 0.02Am _1 , A 2 5 = 0.1 A.m~ 1 , A^Si = 0.05^.m -1 .T -1 . 



this calculation it has been assumed that a single domain of the p— wave component exists 
over the entire length L of the junction. 

However, in real crystals, d+ip domains are expected to occur, whenever the p x compo- 
nent changes sign with respect to the d component, which introduces two other terms that 
were previously cancelled when integrating over the junction. Let us consider n max randomly 
distributed domains along the y-direction. Each domain extends from ^ to J! j^-- Each do- 
main boundary corresponds to a change of sign of the p x component. In the expression for 



the current, Eq. (18), the term proportional to A 2 is replaced by 

-4»e;£s"(-i) ,h - i x 

(*> + SxB) [H (cos(^) - oos(^)) + ft (Bin(^) - sin(^ 



The results of the calculation for multidomains samples are presented on figure 2 for two 
domains in two different configurations (red dotted line and blue dashed line) and 5 domains 
(black solid line). The usual Fraunhofer pattern is replaced by a complicated interference 
pattern. The patterns are not even in field and they do not show a decrease with magnetic 
flux, however the periodicity of a flux quanta is still apparent. (The modification of the 
pattern due to neglected Josephson loop currents introduces a second harmonic but their 
magnitude is expected to be very small.) The observed experimental pattern should vary 
depending on the geometry of the domains, but in every case it should retain a component 
whose magnitude does not decrease with field at large fields. 

When tunneling along the (001) direction of the cuprates, i.e. tunneling perpendicular 
to the Cu-0 planes, the p— order parameter with or without a magnetic field contributes 
no Josephson current. Therefore in this geometry only a conventional Fraunhoffer pattern 
should be observed. 

The situation is more complicated when tunneling along the 100 direction, i.e. along 
the lobes of the d— wave order parameter. A p— wave component would then be necessarily 
accompanied by a current at an angle 7r/4 to the boundary. This is not allowed. A systems 
of domains along the surface can be envisaged which cancel the component of the current 
perpendicular to the surface. We have not investigated the energetics of this possibility, but 
if this is realized, the Josephson current variation with a magnetic field is expected also not 
to be decaying as in a Fraunhoffer pattern, due to the boost of the p— wave by the Meissner 
current. 

We are aware of one set of Josephson experiments in YBa2Cuj,0§+ x with tunneling 
direction perpendicular to the layers as well as along the 110 and 100-direction [Sj. In 
the former the usual Fraunhofer diffraction pattern is indeed observed. In the latter, an 
oscillatory form remains at large $/$o 5 which is consistent with the prediction here and 
qualitatively resembles the simulation on figure 2. The exact pattern is expected to depend 
on the precise geometry of the d + ip domains and therefore is expected to vary only upon 
warming up and cooling down the sample which was indeed observed. However, in these 
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FIG. 2: (Color on line.) Josephson current patterns as a function of magnetic field for a "d + ip" 
superconductor with domains of the p— wave component. The parameters are the same as in 
Fig.(l). Red dotted and blue dashed curves are for two domains of different configuration whose 
borders are given by (-L/2,-L/3,+L/2) and (-L/2,+L/4,+L/2) respectively. The black solid curve 
is for five random domains (-L/2, -0.41L, 0.09L, 0.15L, 0.23L, +L/2). Parameters are A\ = 
0.02A.m- 1 ,A 2 5 = O.lAm -1 , A 2 5x = 0. 05 A.mT 1 .T _1 . 

experiments, it is not known (although probable) whether the thin films were underdoped. 
We suggest experiments in which the quasi-particle tunneling above Tc should be monitored 
simultaneously with the Josephson tunneling below Tc to see if the effects predicted here 
only occur in samples which show the pseudogap in the quasi-particle tunneling. 

This experiment if verified may constitute an additional observation of the symmetry 



8 



breaking in the pseudogap phase of the cuprates. The predictions here also have applications 
in appropriate geometries to other superconductors lacking a center of inversion [7[ . 
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